Rashba spin-orbit coupling and spin relaxation in silicon quantum wells 
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Silicon is a leading candidate material for spin-based devices, and two-dimensional electron gases 
(2DEGs) formed in silicon heterostructures have been proposed for both spin transport and quan- 
tum dot quantum computing applications. The key parameter for these applications is the spin 
relaxation time. Here we apply the theory of D'yakonov and Perel' (DP) to calculate the electron 
spin resonance linewidth of a silicon 2DEG due to structural inversion asymmetry for arbitrary 
static magnetic field direction at low temperatures. We estimate the Rashba spin-orbit coupling 
coefficient in silicon quantum wells and find the Ti and T2 times of the spins from this mechanism 
as a function of momentum scattering time, magnetic field, and device-specific parameters. We ob- 
tain agreement with existing data for the angular dependence of the relaxation times and show that 
the magnitudes are consistent with the DP mechanism. We suggest how to increase the relaxation 
times by appropriate device design. 



^ ; I. INTRODUCTION 

^ . Electron spins in silicon have been proposed as an attractive architecture for spintronics and quantum information 
devices. The inherently low and tunable spin-orbit coupling (SOC) in silicon heterostructures and the possibility of 
ehminating hyperfine coupHngs by isotopic purification bodes well for quantum coherent spin-based qubits and spin 
' transport. Early experiments together with theory have shown that coherence times can be upwards of three orders 

of magnitude longer than in GaAs.P,0,|^,'^ 
- Energy relaxation of locaHzed spin states has attracted theoretical attention 0, 0, S S| and experimental effort 
^2 I 0,01 for decades, and this activity has recently revived in the context of quantum computation. The idea is to store 
quantum information in the spin of a single electron confined in a semiconductor structure, either attached to a donor 
Q ' atom or confined electrostatically in a quantum dot. Spin transport, also of great interest, encodes information in the 
spin states of an ensemble of electrons. In both cases, electron spin resonance (ESR) measurements of spin relaxation 
' provide a key and available measure of spin coherence properties of electrons in silicon quantum wells, though not a 
<^ [ one-to-one correspondence. Our aim in this paper is to explain some existing ESR results for silicon 2DEGs at low 
If^ • temperatures and to make predictions for future experiments. 

The structures that concern us here are layered semiconductor devices of Si and SiGe. The active layer is the 
i quantum well (QW) that confines the electrons in the growth direction. This layer will be assumed to be composed 
of pure, [001] strained silicon. We shall also neglect any roughness or miscut at the the Si/SiGe interfaces. Devices 
made in this way are commonly referred to in the semiconductor industry as MODFETs and are designed to maximize 
mobility. Figure introduces two example structures. 

Extensive theoretical work has been done on spin relaxation in GaAs and other III-V materials. The developments 
"j^ ■ that began with the theory of D'yakonov and Perel' ^10. .IL. .12] are most relevant for our purposes. These authors 
' found that fiuctuating effective magnetic fields due to momentum scattering in the presence of SOC is the dominant 
I ' spin relaxation mechanism in semiconductor 2DEGs for low temperatures. Here, we start with this assumption and 
use a general spin-density matrix approach to calculate the relaxation times of a 2DEG in the presence of a static 
Ch , magnetic field, including explicitly the angular dependence. Understanding the angular dependence of the linewidth 
O ■ is important for comparison with ESR experiments and the extraction of physically relevant parameters such as the 
momentum relaxation time. We also calculate separately in a fc-p formalism the Rashba spin orbit coupling parameter 
in silicon quantum wells. This is a key parameter in our calculation as well as for other spin control considerations, 
both negative and positive. 

In the next section we discuss the origin and magnitude of the SOC in realistic heterostructures. The section 
' following that presents our calculation. Lastly, we compare with experiment and discuss the impHcations for device 
design. 



II. SPIN-ORBIT COUPLING 



The strong macroscopic electric fields inside heterostructure QWs are important for understanding SOC, especially 
in silicon. These fields are also device-specific, so we carry out our calculations on the two representative structures 
in Figure ni Both devices have square QWs, with equal barriers on the top and bottom interface. The first is typical 
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Figure 1: Strong, internal electric fields are common in silicon quantum well devices. Bl: A typical, high-mobiliy SiGe 
heterostructure uses a donor layer to populate a hig h-density 2DEG. The charge separation results in an ~ 10® V/m. 
1)2: A proposed quantum dot quantum computer which utilizes a tunnel-coupled backgate to populate the quantum well 
without the need for a nearby donor layer. Here, i7z > 10^ V/m due to the image potential formed on the backgate. 



of MODFETs and employs a donor layer above the QW in order to populate it. This charge separation produces an 
electric field between the two layers (across the barrier or spacer layer) which can be approximated by 

E, w = -6 X 10^ V/m, (1) 

where ris = 4 x 10^^ m~^ is the density of electrons in the 2DEG for Device 1, e is the charge of an electron, and 
ei are the dielectric constants. We assume that the QW is populated only by donor-layer electrons, leaving an equal 
amount of positive charge behind. The second structure is one that has been proposed for use in a quantum computer 
device. It utilizes a near- lying, tunnel-coupled backgate 2DEG (< 30 nm away) together with Schottky top-gates 
to populate the QW selectively. This situation also results in a strong electric field due to the image potential on the 
back gate. For one qubit, this can be estimated as 

~ -A - 3 X 10^ V/m, (2) 

where d = 20 nm is the distance from the QW to the back gate for Device 2. Schottky top- gates and other device 
parameters can augment or reduce this growth-direction electric field nominally up to the breakdown field of silicon, 
3 X 10^ V/m, or the ionization energy of the electron. Indeed, this field can actually be smaller than that due to 
the top-gates in certain dot configurations. 

The shift of the electron g-factor from its free-electron value 50 — 2.00232 is one measure of SOC in a system. It 
is quite small in bulk silicon and depends on the magnetic field direction in the (elliptical) conduction band minima 
(Agii w -0.003, Ag± w -0.004). However, it is difficult to reliably extract the SOC strength in a 2DEG from Ag. 
Many parameters (e.g., strain, barrier penetration, Ge content {gce — 1-4), non-parabolicity of the band minima) 
infiuence the magnitude and sign of Ag and it may show considerable sample dependence. The non-parabolicity 
effects are especially sensitive to the electron density within the QW and can hide the magnitude of SOC within a 
system. (T^ 

In these silicon heterostructures, SOC is dominated by inversion asymmetry within the device. The spin-orbit (SO) 
Hamiltonian to first order in momentum is given in an arbitrary electrostatic potential V by 

Hso = , \ ^ E^a • (z X p), 

where Gi are the Pauli matrices. Note that the effective magnetic field that acts on the spin is in the plane of the layer. 
In Si heterostructures, the macroscopic fields, which do not average out, are more important than the atomic electric 
fields. In the noncentrosymmetric III-V materials such as GaAs, this is not necessarily the case and the resulting 
Dresselhaus or hulk inversion asymmetry fields are usually dominant. The asymmetry considered here, due either to 
an interface, charge distribution, or external potential, is usually called Rashba or structural inversion asymmetry. 
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The Rashba term comes directly from the SO Hamiltonian if we assume one, dominant symmetry-breaking electric 
field in the structure and average over a momentum state. In a QW, as we have pointed out above, the electric field 
is in the growth (z) direction and thus the z-component of the above dot-product is selected and we obtain 

H]f ^ aipr^ay - pya^) (X E^{a X p)^, (3) 

which is then the Rashba-Bychkov Hamiltonian . [l^ 

Strictly speaking, as de Andrade e Silva et. al. point out ^3|i the conduction-band-edge profile, Ec, and the space 
charge separation (or appHed electrostatic field), E^, contribute separately and sometimes dissimilarly to the SOC. 
For example, the wavefunction discontinuity (band offset) across a material-interface can cause Rashba spin-splitting 
itself. However, in devices of the type considered here, the macroscopic field should be the main contribution. These 
same authors have derived an expression for a in the Kane model for GaAs. We have adapted their work for Si, using 
a 5-parameter 8-band Kane model. This is 8 bands including spin, which means just the lowest conduction band and 
the three highest valence bands. By calculating the breaking of the degeneracy between the spin-up and spin-down 
states of the lowest conduction band, we find 

where we have taken the average of the electric field in the z-direction. Here P = h{X\px\S) /im, = h{Z\pz\S) /im, 
Ad = 0.044 eV is the spin-orbit splitting of the two highest conduction bands, Eyi = 3.1 eV is the direct gap of 
the strained sample, and Ev2 = 7 eV is the gap between the conduction band minimum and the lowest of the three 
valence bands. (These are the 5 parameters mentioned above.) m is the bare electron mass. The matrix elements 
that define P and Pz are to be taken between the cell-periodic functions of the indicated symmetry at the position 
of the conduction-band minimum. Unfortunately, these are not well known in Si, since other bands contribute. We 
may note that P and Pz are examples of momentum matrix elements that don't vary too much in HI-V materials 
and Ge,^3| and are given approximately by 2mP^^^/h^ « 22 eV. With these values we find that for devices of type 1, 

aKilMxlO-^{Ez) m/s. (5) 

Previous Kane models for GaAs involve matrix elements at A: = 0. Our theory is new since it takes into account 
the proper symmetry of silicon with it's minima well away from the zone center. 

Wilamowski et. al. using conduction electron spin resonance (CESR), have measured a « 5.94 m/s (ae/h — 
0.55/v^ X 10-12 gY 

•cm in their units) in a Sio.75Geo.25/Si/Sio.75Geo. 25 QW where the strained-silicon layer was 
roughly 12 — 20 nm. Carrier concentrations were rt^ «:! 4 x lO^^m-^. These numbers correspond to our Device 1 
parameters. Our equations then give, using Eq. 

w 5.1 m/s for {E^) 3 x 10^ V/m. (6) 

Theory compares in order of magnitude and we believe that our estimation has some utility as a guide for device 
design. 

For Device 2, a^^ « 0.25 m/s for (Ez) — 1.5 x 10^ V/m. This device remains to be built. For Device 1, we can 
also predict the zero magnetic field spin-splitting in a silicon 2DEG using |e+ — e-| < 2apF, 



2apF = 2ahJ « 0.75 ^eV, 

where 4 is the degeneracy factor in siHcon (spin+valley). This has not yet been directly measured to our knowledge. 
Taking a Zeeman splitting of g^sB — 0.75//eV with a g-factor of 2, this impHes an internal, in-plane, effective magnetic 
field — the so-called Rashba field — of roughly 62 Gauss, which is the direct result of SOC in the silicon 2DEG of Device 
1. 

Let's consider the relevance of our silicon SOC results to QC and spintronics. We note that the magnitude of the 
Rashba coefficient is much smaller than for GaAs. For the same electric field and 2DEG density as Device 1, a similar 
k ■ p theory for GaAs arrives at aoaAs ~ 230 m/s.|0| But GaAs itself is not a high Rashba HI-V semiconductor and 
is thought to be be dominated by Dresselhaus SOC {0£nAs ~ 1000 m/s).[23| InAs-based heterojuctions, for example, 
may have orders of magnitude higher Rashba values. (2^ This means that SOC effects in siHcon devices will be much 
smaller, including decoherence and gating errors that are SOC-based. 
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Figure 2: ESR linewidths in Gauss as a function of static magnetic field direction (where 9 = is perpendicular to the 2DEG 
plane) for specific values of 2DEG density and Rashba asymmetry. The quantum well is assumed to be completely donor-layer 
populated and as such, a is calculated directly with Eq. |^ as a function of the 2DEG density. 



III. SPIN RELAXATION 



We wish to consider the combined effects of the SO Hamiltonian 



H 



2D 



X <Jy Py (7x ) 



and the scattering Hamiltonian. The scattering may be from phonons or from static disorder. We take the semiclassical 
approach, in which the effect of scattering is to cause transitions at random intervals from one wavepacket centered at 
p with sp = £f to another centered at p* with = ep, where £f is the Fermi energy. This corresponds to a random 
switching in the direction of the effective magnetic field that acts on the spin degree of freedom. This is the D'yakonov- 
Perel' mechanism of spin relaxation. (Tfl The measured quantity in the continuous-wave experiments carried out on 
2DEGs is T2, the transverse relaxation time, while pulsed experiments can also measure Ti, the longitudinal relaxation 
time. For our purposes, a density matrix approach is the natural one, since we will eventually want to perform an 
ensemble average over all possible scattering sequences. Since the physical model of spins in a random time-dependent 
magnetic field is the same as that for relaxation of nuclear spins in liquids, the Redfield technique may be used. 

We outline the calculation only briefiy, since the details are parallel to the discussion in standard texts. [22| The 
2x2 density matrix p allows us to compute the expectation values of the spin by ((Ti) = Tr {(Jip\ ■ For a single 
system described by the Hamiltonian H, we have the equation of motion dp/dt = | [p, H] . In the case of 



(7) 
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Figure 3: ESR linewidth lifetime T2 from Eq. El for constant asymmetry coefficient, a — 1 m/s, as a function of 2DEG 
mobility, fj,, and density, n^. For donor-layer populated quantum wells, divide the times listed by a^: T2{a) = T2{a = 
The magnetic field is assumed to be -B = 0.33 Tesla, perpendicular to the plane of the 2DEG. 



where Hi is small, it is convenient to go to the interaction representation 

pint — ex.p{—iHQt/h)pe:sq){iHot/h), 

and then we get 



dp' 



where 



dt h 



Hl^^t) = exp{iHot/n)Hicxp{~iHot/h). 



This equation can be integrated to give 



and this can be solved interatively, which in second order gives 

i 

h 



= ^[p-(o),i7r*w] 



(8) 

(9) 
(10) 

(11) 

(12) 



■ \ 2 rt 



dt' 



For example, let the steady field be in the z-direction, so that Hq — fkUcaz/^. The fiuctuating field [^i"*(i)] ^^z, — 
J2i=x y^^i^)'^ss" is i'^ the transverse direction. The first order matrix element vanishes and we are left, in second 



order, with 



dt 



J2 P-"(0) [Hr\t')]^„^,„ [iJ™*(i)],,„,,e^(^"-^"')*'e^(^"'-^')* 

S S 

+ J2 [^r*W],,.e'^'~'"^* [iJ™*(t)]^„^,„e'(^"-^"')*Vs"'s'(0) 

s"s'" 
s" s'" 



(13) 
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where we define e'** = e*^'*. Averaging the random field, we find 



where t = t — t' and we define x*(t) = h^{t)h^{t'). We substitute this expression into Eq. Eland do the matrix 
algebra. Eq. El can then be integrated in the limit where t is large. We neglect the oscillating terms, which then 
yield 



= p+-(0) [x^M + xn'^c) + 2x^(0)] 



,h) \ dt 
and 



^hj \ dt 

To obtain the relaxation times we must consider the equation for the spin: 

d(cr^) _ d[Tr{a^p)] _ d 

~ir - — It — -jt^p^+'P-^ 

= 2('-\ [x^M+xn'^c)][p++(0)-p-(0)] 



and so 



Also 



-2(-) [x^M+xn^c)]M (14) 



l/Ti = 2[x^(^,)+X^(^c)]//i' 



which gives 



This gives a relation 



d(cr"^) drr(cr^p) d. 

= [p+_(0) + P+-(0)] r(c..) + x^(^c) + 2x^(0)] /h' 
= - [X"(^c) + xH^c) + 2x^(0)] (^.) //^^ 



1/T2 = [x^(^,)+X^(c.'c) + 2x^(0)] 



I/T2 = I/2T1 + 2x^(0)7^2 



We now wish to specialize to the case of a DP mechanism in a 2DEG. The main point is that the static field may 
be in any direction, while the fiuctuationg field is in the plane. Consider first the special case that the static field is 
along z. Then 

1/T2 = i/2ri = [x^(c.,) + x^(^c)] 

Now consider a general direction, say B along the direction BxX + BzZ — sin9x + cos 9z, so that 9 is the angle to the 
normal. Then the longitudinal fiuctuations x"; which are quadratic in the field, are proportional to sin^ 6* and the 
transverse ones X''"to cos^ 9. Thus 

1/T,{9) - 2 [cos2 0x"(^l) + xH^l)] /h^ (15) 
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while 



(16) 



For the DP mechanism in a 2DEG the random field is constant in magnitude, but random in direction. The statistics 
of this field are Poisson: namely that if the x-component the field at time i = is ft-^, then the chance of it remaining 
at decays as exp(— f/rp). Hence 



In these formulas, Tp is the momentum relaxation time. Note that these formulas assume s-wave scattering. 

The zero-frequency limit of these formulas agrees with the recent results of Burkov and MacDonald ||23] [see their 
Eq.(17)] . They do not agree with the formulas in Wilamowski et. al. [see, e.g., their Eq.(3)] who state that the 
relaxation from the DP mechanism should vanish when 9 = Q. This is not consistent with our results. 

The DP mechanism has the nice feature that it is relatively easy to isolate experimentally. It is strongly anisotropic 
in the direction of the applied field compared to other mechanisms. To illustrate this we plot the ESR linewidths 
as a function of field angle in Fig. El What is most striking is the opposite dependence on angle for the rates 
l/Ti and 1/12, with I/T2 maximized when the field is in the plane of the 2DEG, while 1/Ti is maximized when 
the field is perpendicular to the plane of the 2DEG. Physically, this comes from the fact that the electric field is 
perpendicular, so that the fiuctuations of the effective magnetic field are in the plane. Longitudinal relaxation (Ti) 
is due to fiuctuations perpendicular to the steady field, while transverse relaxation {T2) is due to fiuctuations both 
perpendicular and parahel to the steady field. This mechanism has the characteristic that the change in l/Ti as the 
field is rotated through 90 degrees is always a factor of two. The change in I/T2 is frequency- and lifetime-dependent, 
with the anisotropy increasing as the mobility increases. 

The DP relaxation also has the counter-intuitive inverse dependence of the spin relaxation time on the momentum 
relaxation time: l/ri,2 (x Tp, for small Tp (or zero field), typical for motional narrowing. We plot the dependence of 
T2 on the mobility in Fig. At high mobilities and high frequencies lo » l/r^, we find l/ri^2 oc 1/rp. 



We have calculated the transverse and longitudinal relaxation times of a silicon 2DEG in an arbitrary static 
magnetic field. To test our calculations, we compare them to known ESR data.P,|3,E,l3,l23| We hmit ourselves to 
low temperatures, e_F ~ 10 — 15 K, and realistic material parameters for state-of-the-art heterostructures. 

The most robust prediction of the theory is the anisotropy, particularly that of Ti , which is completely independent 
of ah parameters. The only measurement, in Ref. [3], gives satisfactory agreement for Ti : y^^ll^/y^sj_z _ q qj 
opposed to the prediction 0.5. Furthermore, the anisotropy of I/T2 goes in the opposite direction, as it should. The 
magnitude of this anisotropy is measured to be r^"^/r;^^^ = 12.5 which is about a factor of six larger than the 
theory predicts for the quoted mobility. The relaxation times, as far as can be determined by the range set by the 
uncertainty in silicon band parameters, are in rough agreement with what one gets from estimates assuming that 
this well is a device of type 1. The anisotropy of the 2DEG ESR linewidth is independent of a and is indeed only 
dependent on one free variable: the momentum relaxation time Tp which we assume is directly proportional to the 
mobility. The magnitude of the relaxation time, on the other hand, is set by the Rashba coefficient together with the 
Fermi momentum. 

Table n details results for four other samples as well, also nominally of type 1 for which measurements of T2 have 
been performed. Agreement is good for the sample of Ref. jlill. particularly for the anisotropy of T2. This is a 
donor-layer-populated sample measured with CW-ESR. Ref. [2^, measured via electrically-detected ESR (ED-ESR), 
also is in good agreement with the anisotropy predicted by the theory. Comparison with the sample of Ref. also 
seems to be very good. This is an IBM 2DEG with a density of roughly 4 x 10^^ m~^ and a quantum well thickness 
of 10 nm, fully donor-layer populated. Agreement is considerably less good for the sample in Ref. 0[. Here the 
mobility is not well-known, and the sample is partially populated by illumination. So we are lead to believe that there 
is some difference between the two sets of samples that causes one set to have a larger T2 anisotropy than our theory 
predicts, even for very similar material paramters (density and mobility). Further experimental work along the lines 
of measuring the density, mobility, Ti, and T2 on the same samples is needed. 



h^iO)¥iTp)^{hl)e-'^^ 



2 2 

= -a Ppe 



t/Tp 



and 




IV. DISCUSSION 
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Source 


Linewidth 


Anisotropy 


Ref. [12] 
5-30 K 

M ~ 20 mVV-s 
ris ~ 4 X 10^'' m-^ 

CW-ESR 
donor populated 


Exp.: 

Tj^ll" = 420 ns (0.15 G) 
Ta^-L^ = 140 ns (0.45 G) 
Fred.: 
ra^""" =191 ns 
7f -L^ =60 ns 


yi3||.^yS±Z^ 3 
ySlk^yBXz^ 3.2 

Tf "7rf 0.5 


Ref. [2] 

5 K 
/i ~ 9 mVV-s 
n, ~ 3 X lO^^m"^ 
Pulsed-ESR 
light populated 


Exp.: 

Ta^"^ =3 /iS 
Ta^-^'^ =0.24 MS 
rf =2 MS 
Tf =3 MS 

Pred.: 
Ta^"'' =502 ns 
Ti^' =272 ns 
rf ""^ =251 ns 
=502 ns 


ysil.^ysx.^ 12.5 

Tf "7rf-^'==o.67 

ySlk^yBXz^ 1.8 

rf ii7rf-L^=o.5 


Ref. [1] 

M ~ 10 mVV-s 
~ 3 X 10^5 
CW-ESR 
light /gate populated 


Exp.: 

Ta^"^ =12 MS 
Ta^^^ =500 ns 

Fred.: 
Ta^"^ =480 ns 
Ta^^^ =193 ns 


yS||./yB±.^ 24 
yS||.^yB±Z^ 1.9 

rf "7rf 0.5 


Ref. [24] 

/i ~ 5 m^/V-s 
n« ~ 4 X 10^^ 
CW-ESR 
donor populated 


Exp.: 

Ta^"^ = 50 ns (0.13 G) 
Ta^-L" = 30 ns (0.215 G) 
Pred.: 
Ta^"" =31 ns 
Ti^^"" =20 ns 


ySll.^ySX.^ 1.65 

ySllz^yBX.^ 1.6 

rf "7rf 0.5 


Ref. [2^ 
4.2 K 

/i ~ 9 mVV-s 
Ms ~ 4 X 10^5 m"^ 

ED-ESR 
donor populated 


Exp.: 

Ta^"^ = 105 ns (0.6 G) 
ysx^ = 49 ns (1.3 G) 
Fred.: 
Ta^"^ =212 ns 
Ta^^" =115 ns 


ySllz^yBX.^ 2.2 

ySlk^yBXz^ 1.8 

Tf "7rf 0.5 



Table I: We calculate relaxation times assuming that the QW is completely populated by the donor layer. Accurate analysis is 
made difficult due to lack of precise values for mobility and density, which are often not measured directly (or reported) for the 
specific sample addressed with ESR. The anisotropy does not depend on the Rashba coefficient. Note also that there is some 
disagreement in the literature as to how to convert from linewidth to a relaxation time, we use the equations derived by Foole 
in Ref. |2^, Ta = 27 (^VSg/iB^Hpp) , but others may differ by a factor of up to 2tv. 

The magnitude of the predicted relaxation times is in general well predicted by the theory, but there is a range of 
error. The position of the ionized centers which populate the well is important for the calculation of the electric field 
which is thus hard to characterize. Photoelectrons created with light at the bandgap energy may also have symmetry 
changing effects. This may explain why the Rashba coefficient derived from varying the 2DEG density by light in Ref. 
[Tfll | appears to be independent of density. It is also important to point out that parafiel conductivity (current paths 
through both the 2DEG and the donor-layer for example) is a common problem in today's SiGe quantum wells and 
may eflFect the transport measurements of density and especially mobility making comparison with theory difficult. 

Other mechanisms may become important as we leave the parameter range considered in this paper. Electron- 
electron collisions which do not greatly affect the mobility at low temperatures may start to contribute at higher 
temperatures and mobilities, as they appear to do in GaAs quantum wells. These collisions will also relax the 
spin, but the relation between momentum relaxation and spin relaxation is not expected to be the same as for the 
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elastic collisions considered here. At higher magnetic fields, the cyclotron motion of the electrons is important. In 
the semiclassical picture, when uJcTp > 1, the average value of the momentum perpendicular to the magnetic field 
shrinks, reducing or even ehminating DP spin relaxation, as has been considered for III-V semiconductors in the 
fixed magnetic field case.(l3. [2^ This effect may become important at high mobilities and would be dependent on 
magnetic field angle, increasing the anisotropy predicted here while also increasing the relaxation times. Quantum 
effects may become important in this regime however. The wave vector dependence of the conduction band electron 
g-factor may also lead to relaxation, as has been pointed out recently,!^ but hasn't been considered here. Finally, 
the addition of details related to the presence of two conduction-band valleys may differentiate further the case of Si 
from that of GaAs. Golub and Ivchenko have considered spin relaxation in symmetrical (a=0) SiGe QWs, where 
valley domains (even or odd monolayer regions of the QW) may have influence over spin dynamics. Random spin- 
orbit coupling due to variations in the donor layer charge distribution may also be important in symmetric quantum 
wells. jH Hi 

Long spin relaxation times on the order of hundreds of nanoseconds to microseconds, found in presently available 
SiGe quantum wells, hold great promise for both quantum information processing and spintronics. Our results demon- 
strate that decreasing the reflection asymmetry within the device will appreciably decrease the Rashba coefficient and 
the consequent spin relaxation at low temperatures. This can be achieved by a symmetric doping profile or using an 
external electric field to cancel out the field of the ions. They further show that the anisotropy of the ESR linewidth 
as a function of angle may be a good indicator of 2DEG quality (mobility) independent of transport measurements. 
As higher mobility and more exotic SiGe heterostructures are grown and characterized, new physics may emerge. 
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